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Abstract 



A search for the decay K$ — > 7r°77 has been made using the NA48 detector at the 
CERN SPS. Using data collected in 1999 during a 40-hour run with a high-intensity 
Kg beam, an upper limit for the branching ratio BR (Kg — > 7r°77, z > 0.2) < 
3.3 x 10~ 7 has been obtained at 90% confidence level, where z = ml /rn? K . 



1 Introduction 

The study of the decays K S L — > 7r°77 is useful to test the predictions of the Chi- 
ral Perturbation Theory (%PT) on branching ratios and g 2 -spectra [|TJ. At present the 
branching ratio for the decay Kl — > 7r°77 has been measured precisely [[|, while the 
corresponding Ks decay has not yet been observed. xPT predicts a branching ratio for 
Ks — > 7r°77 of 3.8 x 10~ 8 with the cut-off z = m 2 7 / 'm 2 Ko > 0.2, to avoid the region with the 
pion pole which dominates the total rate fl| . The amplitude can be calculated precisely as 
it is non-vanishing in lowest order Chiral Perturbation Theory and expected to be quite 
insensitive to higher order corrections. With sufficient statistics, the chiral structure of 
weak vertices could be tested experimentally in the high-g 2 region from a comparison with 
the predicted g 2 -spectrum. 

2 Experimental set-up and data taking 

The NA48 detector, originally designed to measure the direct CP-violation param- 
eter e'/ e in neutral kaon decays with simultaneous Kl and Ks beams, was used to collect 
data for studies on rare Ks decays during a test run in 1999 with a high-intensity Ks 
beam. This beam had an intensity about 200 times higher than the one used in the e'/e 
runs. The present analysis is based on data collected during this 40-hour high-intensity 
Ks run. Neutral kaons were produced by a primary beam of 450 GeV/c protons, extracted 
from the SPS accelerator during a 2.4 s spill every 14.4 s and impinging with an angle in 
the vertical plane of 4.2 mrad on a 2 mm diameter, 400 mm long beryllium target, at an 
intensity of ~ 6 x 10 9 protons per pulse. The target was followed by a sweeping magnet 
and a 0.36 cm diameter collimator to define a narrow beam of neutral kaons and hyper- 
ons. The beginning of the fiducial volume was set 500 cm downstream of the centre of the 
target, 1 m before the end of the collimator. The fiducial volume was contained in an 89 
m long, evacuated steel cylinder and followed by the main NA48 detector. The neutral 
particles of the beam were contained in a 156 mm diameter carbon-fibre pipe linked to 
the evacuated cylinder and passing through the detector, in order to avoid interactions 
of neutrons and photons from the beam. A detailed description of the detector and beam 
lines can be found elsewhere ||. 

In the present analysis the following sub-detectors were used: a magnetic spectrom- 
eter composed of four drift chambers, two upstream and two downstream of a dipole 
magnet, for vetoing charged particles; seven anti-counter rings (AKL) of iron and plastic 
scintillator, for vetoing photons outside the acceptance of the main detector; a quasi- 
homogeneous electromagnetic liquid-krypton calorimeter (LKr), to measure energy, po- 
sition and time of electromagnetic showers generated by photons; a sampling hadron 
calorimeter composed of 96 steel and scintillator planes. The electromagnetic calorimeter 
is composed of 13212 readout tower cells of 2 x 2 cm 2 cross section and 27 Xq depth each. 
The energy resolution is 0] 



<t(E) 0.09 0.032 

— - — - ~ 6D 

E ~ E Je 



0.0042, (1) 

with E in GeV. The position and time resolutions for a single photon are better than 
1 mm and 500 ps, respectively, for energies greater than 25 GeV. The overall energy scale 
was determined by a fit of the position of an anti-counter (AKS) used during e'/e runs. 

Events satisfying the trigger conditions for 77 or 3n° events were selected, since no 
dedicated trigger for 7r°77 (and 27r°) decays was set up during the special high-intensity 
Ks run. The trigger decisions were based on quantities reconstructed from orthogonal 
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zv < 9 m; npx,npy=l,2 
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npx > 4 or npy > 4 
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~ 1.1 x 10 6 



Table V. Conditions and downscaling factors D for 77 and 3tt° triggers, for different 
parts of the run. zv= z^ftex — 6 m is the position of the decay vertex starting from 
the beginning of the decay region, npx and npy the number of energy peaks in x and y 
projection respectively, in a time window of 20 ns. In the rightmost column the number 
of 27T° events after all selection cuts satisfying 77 or 3tc° trigger conditions is given. 



projections of the energy deposit in the electromagnetic calorimeter |J. The trigger for 
both 77 and 3ti° decays required: 

1. the total deposited energy E tot = E LKr + E H ad being larger than 50 GeV; 

2. the centre of gravity of the event, computed from the first moments Mi >x and Mi tV 
of the energy in each projection, 



co.g tn9 = -, (2) 



being less than 15 cm from the beam axis; 

the decay vertex of the K°, computed from the second moments M^ x and M2 >y , 



„ E LKr (M 2 , x + M 2tV ) - (Ml x + Ml 

trig _ _V v ZJ_ (o\ 

^vertex ^LKr j \°) 

m K 

being less than 15 m downstream from the end of the collimator. This corresponds 
to about 2.5 Ks lifetimes, since the K$ have a mean energy of 110 GeV. In equation 
0, ZLKr is the z coordinate of the electromagnetic calorimeter with respect to the 
target and mx is the kaon mass. 
The 37T° trigger required also at least 4 energy peaks in each projection. The cuts on the 
vertex for the 77 trigger and on the maximum number of energy peaks for the 3ir° trigger 
were varied during the run, while the downscaling factors were adjusted according to the 
obtained rate. The different trigger conditions and downscaling factors are summarised 
in Table |1| 

The trigger efficiency for Ks — > 2tt° decays was computed using a minimum bias 
trigger which did not involve the electromagnetic calorimeter: for the first part of the run 
(conditions a-d in Table |) it was (98.52 ± 0.06)%, while it decreased to (22.0 ± 0.2)% for 
condition e. About one third of the statistics used for this analysis was collected under 
this last condition. The trigger efficiency for Ks — > 7r°77 events was assumed to be the 
same as for 27r° events. Possible differences were studied from simulation and considered 
as systematic uncertainties. 



3 Event selection 

To select Ks — > 7r°77 candidates, events with at least four electromagnetic clusters 
were considered. Each cluster was required to satisfy the following conditions: 
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1. the energy had to be greater than 3 GeV and less than 100 GeV; 

2. each cluster had to be well inside the calorimeter acceptance, more than 15 cm and 
less than 113 cm from the axis of the beam pipe; 

3. the distance from any dead cell of the calorimeter had to be greater than 2 cm; 

4. the distance from any other cluster in the event with an energy greater than 1.5 GeV 
had to be greater than 10 cm. 

For each event, all combinations of four clusters passing the pre-selection described above 
were considered. For each combination, the following selection criteria were applied: 

1. the time of each cluster (t c i) had to be within ±3ns of the event time (t aV 4) computed 
as the average of t c i weighted over the energy; 

2. no other cluster (not belonging to the combination considered) with an energy 
greater than 1.5 GeV had to be found within ±5ns of the event time, to reduce the 
effect of accidental background; 

3. the sum of the four cluster energies had to be greater than 70 GeV and less than 
170 GeV; 

4. in order to distinguish 7r°77 decays from 2ti decays, a x 2 variable was defined f5|: 



X2n° 



^m 77 + m 77 ) /2 - m n o 



+ 



(4) 



where m 77 and m 77 are the invariant masses of the two 77 pairs and a± the res- 
olutions of (m 77 ± rn'^ /2 measured from the data. For each combination of four 
clusters the 77 pairs giving the lowest x 2 were considered and, if more than one 
combination of four clusters satisfied the first three requirements, the one with the 
lowest x 2 was selected. 
Further requirements were then applied to each event: 

1. to define a good event, cuts in the centre of gravity and decay vertex position of the 
kaon were used. The centre of gravity 



v '(EiEiXi) 2 + (EiEiVi) 2 
c -°-9- = , (5) 

with Ei, Xi, yi respectively the energy, x and y coordinates in the LKr for the z-th 
cluster, had to be less than 7 cm; 
2. the decay vertex 
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Zvertex = Z LKr EfEj ■ (Xj - Xj) + (jji - JJj) , (6) 



i,j<i 



was required to be within the first 9 m of the decay volume, in order to reject 
background from Kl — > 37r° events with missing or overlapped photons; 

3. events with signals in the drift chambers within ±10 ns of the event time or with 
overflows in chamber 1, 2 or 4 within ±250 ns were rejected. The overflow condition is 
generated in the drift chamber readout whenever more than seven hits in a plane are 
detected within 100 ns, in order to avoid recording events with high hit multiplicity; 

4. in order to reject events with overlapping showers, an energy- dependent cluster- 
width cut was applied, using a parametrisation taken from K$ — > 27r° decays; 
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5. events with signals in any of the AKL rings within ±3 ns of the event time were 
rejected; 

6. if the energy deposited in the hadron calorimeter was greater than 3 GeV, events 
were rejected if any HAC signal was found in a time window of ±15 ns around the 
event time. 

4 Background rejection 

Background to K$ — > 7r°77 may come from misreconstructed Kg — > 2tt° with four 
clusters in the acceptance; from Kg — > 27r° with three clusters in the acceptance and one 
random cluster accidentally in time; from K$ — > tt ^ with one lost particle and one 
misidentified charged particle; from Ki — > 37r° with two lost or overlapped photons; and 
from Ki —>■ 7r°77. 

4.1 Background from misreconstructed K$ — > 2tt° 

To define a good 7r°77 event, the 77 pair with invariant mass closest to the nominal 
7T° mass was required to be within ±3 MeV/c 2 of the nominal ir° mass, corresponding to 
about 3 standard deviations of the resolution. K$ — > 2tt° events with misreconstructed 
cluster energy due to non-Gaussian tails in the energy resolution were simulated and found 
to be concentrated in a region with x\^ below 300. Therefore a xL° greater than 300 
was required for the 7r°77 candidates. In order to compare data and simulation a control 
region with 50 < x\^ < 300 was defined and the distribution in xln * s shown in Fig. [I] 
left. The cut on x^o discards about 8% of signal events (Fig. |l| right). The remaining 77 
pair had to satisfy the requirement z = vr? 1 /'m 2 KQ > 0.2, in order to reject background 
from 2-7T events and to be able to compare the result with the prediction made by xPT, 
which is given with the cut-off z > 0.2. After all cuts, 2 Ks — > 7r°77 candidates were 
selected, i 

From Monte Carlo simulation the background from Kg — > 2ir° decays with four 
clusters in the acceptance is expected to be less than 0.5 events at 90% C.L., which 
corresponds to a rejection factor of at least 4 x 10 7 . 

4.2 Background from Kg — > 27r° with three clusters in the acceptance and 
one random cluster in time 

Using simulated Kg — > 2n° events and events in the data selected by a random 
trigger, the background from Kg — > 2ir° decays with three clusters in the acceptance and 
one random cluster accidentally in time was estimated to be 3.0 ± 0.1 s t a t events. Random 
events were assumed to have a flat distribution in time with respect to the average time 
of the 3 clusters of the 2ti° event (t av3 ) and were overlayed to the simulated K s — ► 2n° 
with three clusters in the acceptance to reconstruct a K$ — > 7r°77 candidate. One of the 
2 selected Ks — > 7r°77 candidates in the data contains one cluster for which the time 
is compatible with the distribution of the random cluster in background simulation. For 
the other Ks — * 7r°77 candidate, the farthest in-time cluster from the event time has the 
same probability to be a random cluster or to be part of an event with four photons in 
the acceptance of the calorimeter (Fig. |2|). 

4.3 Background from K s — > tt ^ and K L — > 3n° 

Applying all the selection criteria described above to simulated Ks — * Tt°7[° D and 
Kl — > 3tt° events, no candidate remains, corresponding to a rejection factor of at least 
2.5 x 10 6 for K s — > n ir° D and of at least 7.5 x 10 5 for K L — > 3n . The background from 
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Ks — > 7r°7r% is expected to be less than 0.2 events, while the background from Kl — > 37r° 
is expected to be less than 0.1 events. 

4.4 Background from Kl — > 7r°77 

A non reducible source of background is given by the decay Kl — > 7r°77, which 
cannot be distinguished from the signal. Using the branching fraction measured by the 
NA48 experiment 0, the background from this channel was computed to be about 0.1 
events. 

5 Normalisation and acceptance evaluation 

The decay Ks — > 2tt° was chosen as the normalisation channel. Since the two decays 
differ only in the invariant mass of one 77 pair and were collected by the same trigger, 
systematic uncertainties are expected to cancel in the ratio. The same selection criteria 
were applied, with the only exceptions of x^-o, w hich had to be less than 27, and the 
invariant masses of the 77 pairs, which had to be within ±3 MeV/c 2 of the nominal 7r° 
mass, considering the pairing which gives the lowest x 2 - After all cuts, 3.3 x 10 6 K s -> 2tt° 
decays were selected, with an acceptance of 17.3% computed from a Monte Carlo sample, 
which included non-Gaussian tails in the energy resolution of the calorimeter, but had no 
accidentals and no simulation of the trigger. The comparison between data and Monte 
Carlo simulation in the xL° anc ^ z vertex distributions for selected K s — > 2n° events is 
shown in Fig. [| At high values of z vertex a slight inefficiency appears in the data due to 
the decay vertex cut at trigger level. Correcting for this inefficiency, data and simulation 
are in better agreement at high z veriex . The correction is, however, not applied to the data 
in the analysis, since the trigger inefficiency is assumed to be the same for signal and 
normalisation channel. 

The total number of Ks with an energy between 70 and 170 GeV and decaying in the 
first 8 m of the decay region was measured from Ks — > 27r° events to be (4.1 ± 0.2) x 10 s , 
taking into account the losses of statistics due to trigger downscaling (~ 27%) and to 
rejection of events with accidental in-time hits in the drift chambers or in the AKL 
(~ 55%). This flux was also computed, as a cross check, using the decay K L — > 3n°, 
which however is statistically less significant, and found to be (3.9 ± 0.4) x 10 8 . The two 
results are compatible within the error, which is due to the statistics used to compute 
acceptance and trigger efficiency and to the uncertainty on the branching ratios. The total 
number of Ks is not used in the limit computation, since the same trigger efficiency is 
assumed for signal and normalisation channel. 

For the signal an acceptance of 15.6% was determined from simulation where the 
amplitude predicted in xPT jlj was used. The acceptance as a function of z is shown in 
Fig. £|. The cut on z > 0.2 is expected to introduce a negligible systematic uncertainty, 
since the acceptance for the signal is nearly flat in that region of z. 

6 Systematic uncertainties 

The following effects were studied: 
1. the background estimation was checked with an independent method based only on 
data. The background from Ks — > 2tt° decays with three clusters in the acceptance 
and one random cluster accidentally in time was assumed to have a flat distribution 
in the region 3 ns < \t c i — t av4 \ < 6 ns for the farthest in-time cluster and to be the 
only class of background to be found in that region after applying most of the selec- 
tion criteria. The expected background was then extrapolated from the sidebands 
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to the signal region and compared with the expectation given by the simulation, 
applying the same selection. The two results are compatible within the statistical 
error and the relative difference of 24% is computed as systematic uncertainty. 

2. the different cut on xln between signal and normalisation, together with a slightly 
different resolution of the best reconstructed tt° mass between data and simulation, 
produced a systematic uncertainty of 0.6% on the ratio of acceptances when a 
window of ±3 MeV/c 2 is used for the cut on the 77 invariant masses; 

3. from simulation the relative difference in the trigger efficiency between K$ — > 7r°77 
and K s — > 2n° was computed to be about 1.7%; 

4. not including the effect of non-Gaussian tails in the simulation produced an increase 
in the single acceptances of about 5%, while the variation in the ratio was 1.5%; 

5. the uncertainty of 0.9% on the Ks — > 2ir° branching ratio gives a systematic uncer- 
tainty of the same amount; 

6. the statistical error on the simulated sample is about 0.5%; 

7. varying the selection cuts had a negligible effect on the result. 

The total systematic uncertainty was computed as the sum in quadrature of the 
individual uncertainties and estimated to be 24.1%. The effect on the result is an increase 
of the BR limit of about 6% 0. 

7 Result and conclusions 

Applying all the selection criteria, 2 Ks — > 7r°77 candidates were found and 
3.1 ± O.lstat ± 0.7 syst background events were expected, mainly from K s — > 27r° decays 
with three clusters in the acceptance and one random cluster accidentally in time. The 
acceptance was computed to be a^o^ = 15.6% for the signal and a 27T o = 17.3% for the 
normalisation channel, in which iV 27r o = 3.3 x 10 6 events were selected. The same trigger 
efficiency was used for the two decays. Using N w o^ < 3.14, which corresponds to a confi- 
dence level of 90% when 2 events are seen and 3.1 background events are expected with a 
systematic uncertainty of 24% J7L M, an upper limit for the branching ratio of the decay 
K s — > 7r°77 was calculated: 

BR(Ks^ 7r°77,z>0.2) = • ■ BR (k s -> 2ti ) < 3.3 x 10~ 7 at 90% C.L, 

with z = 

K 

This is the first upper limit on this branching ratio to be published. It lies about 
one order of magnitude higher than the branching ratio predicted by Chiral Perturbation 
Theory. 
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Figure 1: Left: xln distribution for data (squares) and background simulation (solid lines) 
in xlir control region, before the cut on z. Distributions are normalised to the Kg flux. 
Right: xL° distribution for data (squares), background simulation for misreconstructed 
Ks — > 2tt° (filled histogram) and Kg — > 7r°77 simulation (solid line). For the Kg — * 7r°77 
simulation the normalisation is arbitrary. 
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Figure 2: \t c i — t av ^\ for the farthest in-time cluster from the event time in 2tt° decays 
in the data (left) and for the random cluster in background simulation (right). For the 
random cluster a flat distribution in \t c i — t av3 \ was assumed. The full square and triangle 
represent max (\t c i — t av4 ^\) for each of the two selected events. 
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Figure 3: xI-k (left) and z vertex (right) distributions of the selected K$ —> 2ti events for 
data (squares) and simulation (solid line). In the left plot the complete selection has been 
applied except the xl^ cu ^ an d the distributions are normalised to the total number of 
events. In the right plot the complete selection has been applied except the z veriex cut 
and the distributions are normalised to the bin with the maximum number of events. The 
data distribution corrected for the trigger inefficiency (triangles) shows a better agreement 
with the simulation at high values of z vertex . 
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Figure 4: Acceptance as a function of z — m? /m 2 K0 for K$ —>■ 7r°77 simulated events 
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